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Abstract. We investigate stability of replica symmetry breaking solutions in 
generalized p-spin models. It is shown that the kind of the transition to the 
one-step replica symmetry breaking state depends not only on the presence or 
absence of the reflection symmetry of the generalized "spin" -operators (j but on 
the number of interacting operators and their individual characteristics. 
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1. Introduction 

In the study of spin glasses the central role plays the Sherrington-Kirkpatrick (SK) 
model Q . It appeared as an attempt to describe unordered freezing of spins in dilute 
magnetic systems with disorder and frustration. This problem was soon solved at 
the mean-field level. It was demonstrated that the glassy phase in the SK-model is 
characterized by the the full replica symmetry breaking Later the "p"-spin 

model was introduced as a natural generalization of the SK-model; p spins interact 
in this model at each point of the lattice. It was shown in Ref. Q that in this case 
of p-spin interactions when {p > 3) the full replica symmetry breaking (FRSB) does 
not appear at the glass transition temperature Tc and the 1-step replica symmetry 
breaking (IRSB) solution is stable in contrast to the SK model with 2-spin interaction. 
In addition it was discovered that the glass order parameter in the p-spin model is 
discontinuous at Tc contrary to the SK-model. 

The use of p-spin glass models as models for understanding structural glasses was 
pioneered in Ref. 0] . It was shown that these models qualitatively describe many 
aspects of the glass transitions in liquids, e.g. two critical temperatures. The number 
of metastable states in these models is similar to that obtained in numerical modeling 
in liquids. The structure of the dynamical equations for the correlation functions of 
supercooled liquids in mode-coupling theory and for p-spin model are similar [|, |. 

Till now p-spin glass model remains to be in the focus of intense investigations 
(see, e.g., the recent papers Ref. |ll|, ^ |l3[ |l^) since it is a good starting 
point to understand the physics of real glasses. Some aspects are still far from being 
completely understood. We consider a generalization of the p-spin model of Ising 
spins where arbitrary diagonal operators U stand instead of Ising spins p^ . 
The operators have different meanings depending on the problem under study. For 
example, Ising spin should be replaced with the molecule multipolc moment if freezing 
of the orientational order is the target of the investigation |l^ . 

An important development of the p-spin model is the study of the quadrupole 
system with J — 2 with multiparticle interaction, p — 3 [|9|. For the full molecule 
momentum, J — 2, and = {0, ±1,±2}, the operator of the axial quadrupolar 
moment of a molecule has the form, U = ^{3J^z — 6). Using the operator U instead of 
the spin in the p-spin model it is possible to describe the orientational glass phase 
observed for high pressures in solid molecular ortho-D2 and para-iJ2- Moreover, 
it is possible to observe orientational transitions in such systems, which consist of 
initially spherically symmetric molecules with J = 0, because the probability of the 
transition, ( J = 0) — > ( J = 2), increases rapidly with the pressure. The computational 
results show that the glass state and the long-range orientational order coexist, which 
agrees with experiments po| . This model well describes high pressures since there the 
interactions of more than two particles play an important role. 

We consider now a generalization of the p-spin model in the following way: 
the role of spins play the operators satisfying the condition of reflection symmetry, 
any integer k. Reflection symmetry of the operators U leads to 
vanishing of a number of terms in the free energy, so that the replica symmetric 
(RS) solution for the order parameters is zero at high temperature. As the result 
the behavior of IRSB solution for the order parameters is like in the ordinary p-spin 
model of Ising spins with (p > 3) 

If operators U do not have the reflection symmetry, Tr J/^^'^+i) ^ 0, then 
the glass freezing scenario is absolutely different from the Ising p-spin case. The 
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characteristic properties of systera develop themselves already in the replica symmetry 
(RS) approximation. The nonlinear integral equation for the RS-glass order 
parameter simply has no trivial solutions at any temperature because the integrand is 
nonsymmetric due to the cubic terms in the free-energy expansion [ p^ , pT[ . There is a 
smooth increase in the order parameters (both glass and regular) as the temperature 
decreases. Thats why IRSB order parameters appear continually at the bifurcation 
point To- [In some sense IRSB solution behaves here like in ordinary p-spin model 
with spins but in the external field p2[.] 

We found mi analytically in the branching point and expressed the result through 
RS-order parameters. If mi < 1 in the branching point then IRSB solution has 
physical meaning near Tq. The quadruple glass with J — 1 and Jz — {0, ±1} (three 
particle interaction) is the simplest example of the system without the reflection 
symmetry. In this case U = 3J^z — 2 is the quadrupolar moment of the molecule 
(see Fig. 1 in Ref. Q). Then IRSB solution appears to be stable and it branches 
continuously at the bifurcation point Tq = TmsB and smoothly on cooling. 

If in the branching point mi > 1 then IRSB solution does not have physical 
meaning in the vicinity of Tq. We suggest below an illustrative example of this 
conjecture. At other realizations of U operators with Tr [/(^fc+i) ^ q the transition 
from RS to IRSB does not take place at the bifurcation point To where IRSB-solution 
appears (see Fig. p. Formally IRSB parameter mi > 1 at To, however only mi < 1 
have physical sense. While mi > 1 then IRSB solution formally exists but it is 
unstable. The temperature TmsB where stable IRSB appears coincides with the point 
where mi = 1. At this point, F^s — i^iRSB- When the temperature is decreased mi 
becomes smaller than one and the IRSB solution leads to a larger (preferable) free 
energy than the RS solution. 

Below we investigate the crossover from continuous to jumpwise behavior of the 
glass order parameters in generalized p-spin models using the bifurcation theory and 
analyze the stability. We show that in general analytical progress can be made in the 
bifurcation region. We expand the boundaries of the "ordinary" p-spin model and 
consider IRSB solutions also for the pairwise interactions. 



2. Generalized p-spin model 

2.1. Main equations 



The Hamiltonian of the p-spin model in general looks like: 

it<i2...<ip 



H — — ^ Ji,...i„Ui,Ui.,...Ui. 



(1) 



where U now is arbitrary diagonal operator with Tr [/ = 0, A'^ is the number of sites 
on the lattice, i = 1, 2, ...A^, and p is the finite integer giving the number of interacting 
particles. The coupling strengths are independent random variables with a Gaussian 
distribution 



P{Jn...ip) = r-i-r exp 



p\ttJ 



p!j2 



(2) 



Using replica approach we can write in general the free energy averaged over 
disorder: 



1 



n^O n 



4 



{F)j/NkT = lirn^ £ max { - > ]{w°'Y + ^ 
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(3) 



(4) 



The extremum in Eq. (||) should be taken over the physical order parameters 
and over the corresponding Lagrange multipliers, X"^ and fi". So the saddle point 
conditions give the glass order parameter g"^ 



Tr 


tf^tjl^ exp ( 
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Tr 


exp (ej 





the auxiliary order parameter w°', the regular order parameter x": 



Tr 


(C/")2exp 






Tr 


exp 









Tr 


J7" exp (ei 
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Tr 
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and the parameters 



(5) 



(6) 



(7) 



To proceed with replica symmetry breaking procedure it is more convenient to 
rewrite Eq. in the form: 



1 f 

{F)j/NkT = lim - max { (p ~ l)-rY' (w"? 

k a 



b - l)lT E - lnTr{^.} exp ^ I . 



a>/3 



where 



(8) 



(9) 



Using the standard procedure (see, e.g., Ref. we perform the first stage of 
the replica symmetry breaking (IRSB) (n replicas are divided into n/mi groups with 
mi replicas in each) and obtain the expression for the free energy. Order parameters 
are denoted by q"^ = ri if a and fS are from different groups and q"^ = ri + vi if a 
and P belong to the same group. So 



i^msB = ~NkT { mit^p - 1)^ + (1 - mi)(p - l)t^ 



1 



/ dz'^ln / ds 

mi 



Trexp 



(10) 



Here 



zt 



pri 



(p-i) 



■U + st\ 



lp[{ri + vi) 



(p-i) 



ri(P-i)] 



U+ 



^, pK(^-i)-(ri+t;i)(P-^)] ^^, 



(11) 
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and 



dz / z 



/27r 



(12) 



The extremum conditions for i^msB yield equations for the glass order parameters ri 
and vi, the additional order parameter wi, the regular order parameter xi and the 
parameter mi : 



Ti = dz 



Jds^ 


Tr exp flmsB 


(mi-l) 


TtU exp 6*^33 


JdsO 


Trexp e^jiss 


mi 



vi — dz 



Jds^ 



Tr exp 6*15 



fmi-2) r 



TrJ/exp 



/rfsG Tr exp Pins 



d2« 



Tr exp cfipi 



(mi-l) r 



Tr C/ exp 6Iirse 



JdsG Tr exp c/iRs 



wi = dz 



xi = dz 



Jds^^ 


Tr exp fliRSB 


(mi-l) 


TrO-^exp^iHSB 


Jds^ 


JdsG h 
Tr exp 6*^33 


r exp 6'ip 

(mi-l) 


- mi 

SB 

Tr U exp ^irsb 



Tr exp c'lH 



and 



mi 



dz"- 



Jds 



Tr exp 6'iB 



In 



Tr exp Pip 



Tr exp 



J ds^ Trexp 6*1^33 



(13) 



(14) 



(15) 



(16) 



(17) 



If operators U do not have the reflection symmetry then the nontrivial solution for 
the regular order parameter x appears in spite of the fact that x is absent in 9 (since 
Jo = (^...J=Oin @) @0. 

The heat capacity can be expressed through the glass order parameters: 



C'v(iRSB) d 



kN d{l/t) 



mir\ + (1 — TOi)(ri + viY — wi^ 



(18) 



The corresponding expressions for the RS approximation can be easily obtained 
from the preceding formulas (p^-(|l5|) by setting t;i = 0. For glass order parameter 
we have: 
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Here 



= zt 



.(p-i) 



■ u + t 



(20) 



2.2. Stability of the mean-field solutions 

The stability of the saddle point solution can be tested by the investigation of the 
gaussian fluctuation contribution to the free energy near this solution. The mean field 
(saddle point) solution is stable while all the eigen modes of the fluctuation propagator 
are positive. The most important mode is the so-called replicone mode ||, |l^ since 
only its sign is usually very sensitive to the replica symmetry breaking degree and 
to the temperature. For example, the replica symmetric solution is stable unless the 
corresponding replicon mode energy A(i,s)„pi > 0. The RS-solution can break at the 
temperature Tq determined by the equation A(RS)ropi — 0, where 



A(RS)ropi - 1 - r / az 



Tre^Rs 



Tr C/e^i^s 
Tr e^Rs 



(21) 



The equation A(rts)^^pi = can be obtained as the branching condition for Eq. (|lj), 
i.e., as the condition that a small solution with IRSB can appear. 

We see that the equation for the glass order parameter, Eq.(|9|), contains 
'Pr[t/e«Rs]/Tre^"^ If operators U have zero trace for aU odd powers: TrC/(2'=+i) = 
for all integer k, then Eq. (|l^) always has trivial solution = 0. 

In fact, since Eq. ( p^ can have other positive solutions, the stability condition 
is very useful to get the physical one. The stable solution for T ^ is the trivial 
solution that bifurcates. Then equation A(p,s)„pi(i) = is solvable only for p = 2 at 
finite t ~ to. Only in this case the nontrivial IRSB solution smoothly branches for 
To ^ . For p > 2 and Tr C/(2fc+i) = 0, we obtain Tq = 0. So IRSB solution at 
TiRSB 7^ appears discontinuously at mi = 1. 

On the other hand, in the case Tr tj^^'''^^'> ^ 0, the high-temperature expansion 
of the equation for the order parameter does not give trivial solution. Using the 
condition, g^s ^ 0, we can find that Tq 7^ 0. Then the solutions with the unbroken 
symmetry may appear continuously not only for p = 2 but also at any integer p. For 
example, the stable continuous IRSB solution exists for p = 3 when tj = 3Jz^ — 2 is 
the axial quadrupole moment in the subspace J — 1 with Jz — 0,±1 (see Fig. 1 in 
Ref. [^). Below we prove these statements. 

In order to get in general IRSB solution near the bifurcation point Tq where it 
is close to the RS-solution we expand the free energy (|l)-(^ up to the third order, 
assuming that the deviations 8q^^ from and p from are small. 

We use the notation Ai^ for the difference of the free energy i^(iRSB) from its 
replica symmetry part F^^^^^^ . So, 

= J 2 9rs W\^-[r- (mi - l)vi] - vi mi(l - mi)| - 

4 , 

r - (mi - l)vif - t^ IC [r - [mi - l)vi]^ + 

D [r — (mi — 1)^1] wi^mi(mi — 1) — i?3Ui'^mi^(mi — 1)-|- 

BWmi{mi - l)(2mi - 1)} + *(p) + ... (22) 



2 
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where t = tg + At, ri = q„s + "i", Winsn = + p and the expressions for W , 
L, C, D, B^, B4 and ^E* are some combinations of operators averaged over the RS- 
sohition (see Appendices A and B). For example, the coefficient L enters AF hke 

hm„^oiE.,^./'5<?"'^<5-?"^. 

First we consider the case 
conditions for the free energy ( p^ and the inequahty L\t=to ^ 0, we obtain the 
branching condition 

r ~ {mi ~ l)vi ^ + o{Atf . (23) 

This condition states that there is no hnear term for the glass order parameters. 
There is no other linear term because [l — t'^W] \t=to = ARSrcpi|t=to = at the branch 
point. Since the coefficients Ak \ta in the expression for 4'(p) (Appendix B) are in 
general nonzero then from the extremum condition we get, p [r — {mi — l)vi] and 
*(p) = + o(At)''. Finally, we get: 

2mi(l - mi)rAt = 3tomi(l - mi) [-B4 + mi{-Bi + 2B4)] vi , (24) 

(2mi - l)rAt tl {(2mi - 1) [-B4 + mi(-B3 + 2B4)] + 

mi(TOi -l)(-B3 + 2B4)}wi, (25) 

where F is given explicitly in Appendix B. 

Here B3, B4 and F are taken at T = Tq. Then we find from (||) and (|25|) (the 
cases mi — and mi — 1 should be investigated separately, see, e.g., Ref. ||25|): 

mi = B4/S3. (26) 

at the branch point Tq where IRSB-solution appears and 

ri = qns + {mi - l)vi, vi ^ At] (27) 

in the neighborhood of Tq. 

The coefficient of proportionality ( p7| ) depends only on RS-solution at Tq: 

VI ^ ^(^-^^ , , |l + + Mp^,(.-)t1 At, (28) 

2 6-84(1 - mi)to5 I 2 g^s 4 tos j , i ; 

where 
T = 



(P-I) , ^0(P-l) .,(p-2) , \ ^ f ip-l) to(p-l) (p_2) , 



(29) 



where q^s — ^^rs and w'^s = ^ifRs- Expressions ii'i and if 2 are quite long and they 
are written in Appendix B. 

IRSB solution appears smoothly in most cases from RS solution. However if the 
order parameter in the brunching point takes the unphysical values mi > 1 then at 
temperature where mi = 1 and physical (stable) IRSB solution appears it is drastically 
different from the RS solution and as the result the jump singularity appears in the 
heat capacity (e.g. Ref. |l^). The free energy does not have the discontinuity in this 
case (see equation (10)). 

As an example we consider the case of quadrupolar glass U = Q — [SJz^ — 6] 
p9| in the subspace J = 2 in detail and write the explicit solutions of Eqs. (|T3|)(p7p 
for the number of interacting particles p ~ 3 (see Fig. |l|). The glass state and the 
state with the long-range orientational order coexist. IRSB solution appears smoothly 
from RS solution. But in this case, the condition A(RS),„pi = does not determine the 
physical solution in the neighborhood of the branch point Tq, namely, a transition to 
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the nonphysical branch of the free energy takes place. In fact, the transition from 
the RS to the IRSB solution occurs jumpwise at the point Tirsb > Tq determined 
by the condition mi = 1. At this point, Fj^g = ^\rsb- Replica-symmetric solution is 
stable above T^^sb ■ When the temperature is decreased mi becomes smaller than one 
and the corresponding physical IRSB solution corresponds to larger (preferable) free 
energy than the RS-solution. An exceptionally important property of the model is 
that there exists a domain of stability where the IRSB-solution remains stable under 
the further RSBs (see below). 

Now we shall consider the case TrC/^^'^+^) = and investigate IRSB-solution 
near the branching point. Let us return to the expression ( ^2|) for free energy at p = 2 
since only in this case the nontrivial IRSB solution smoothly branch at finite t. Note 
that in the case of zero RS-solution for the glass order parameter the expansion does 
not contain the terms where some indices occur only once. In the case of reflection 
symmetry operators there is no terms where some indices occur odd number of times. 
L\t=to = Oj (see. Appendix B), the branching condition ( p3| ) for IRSB fails. Moreover 



at the branching point we have: B4 



Bo 



B' 



A. 



A12 



C = 2i?3 D = — 3B3. From the extremum condition we find ^(p) = 
mi — 0; 



Al4 

4 



ri + Vi 



o{At) 



At. 



= 0; 
and 

(30) 
(31) 



{U?UiUi)to' 

It have been shown that the Parisi FRSB scheme can be used not only for the 
SK model but also for any model with pair interactions and with reflection symmetry 
j'j-jji^k+i) _ example, such as spin glasses with arbitrary spin. 

We break the RS once more and obtain the corresponding expressions for the free 
energy and the order parameters. The bifurcation condition A(ijisB),„pi = determining 
the temperature T ~ T2 follows from the condition that a nontrivial small solution for 
the 2RSB glass order parameter appears as V2 0. We have: 



A 



(iRSB)ropl 



l~t 



,p(p- l)(ri -f wi)(P-2) 



Jds^ 



Tr exp 



Tr[a^ cxp(eiRSB)] 
Tr[cxp(eiRSB)] 



Tr[Lrcxp(eiRSB)] 



n 2 



Ti-[cxp(9irsb)] 



jdso 



Tr exp 



^^IRSB^ 



(32) 



Note that Eq. (32) depends only on IRSB-solution. The expression for A 



(lRSB}rcpl 







always has the solution for wi = 0, which determines the point Tg ^.nd coincides with 
the solution of Eq. ( ^ ) \(T^s)rapi — (see Fig. |^). 

Using the expressions for the glass order parameters obtained above near Tq we 
can show that ^A(iRSB)rcpi Wo— ss for the two-body interaction in the presence of 
reflection symmetry and for the caseTr [7(2*^+1) ^ Q and arbitrary p. 

In addition to the point Tq, one more bifurccition point A(^iptg3^j.(.pi — (see, e.g. 
Fig. H) may exist as vi ^ 0, and the 2RSB solution can appear at this point. At the 
point T2 a transition to FRSB-state or to a stable 2RSB-state may take place. 

We now present several results that hold in the two cases above. They concern 
the general form of the expressions determining the stability of IRSB solutions. Now 
we will investigate the stability of the IRSB-solution to small perturbations having 
2RSB symmetry for any value of the inverse temperature t. Group of replicas of mi 
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Figure 1. (Color online) Order parameters evolution with the temperature for 
3-quadrupole model in the subspace J = 2. For simplicity only the physical 
solution for xirsb is shown here. There are four characteristic temperatures in 
the model: T2, To, Tirsb and T* , where Ajirsbj = 0, Aj^g) = 0, miRSB = 1 and 
the branches of IRSB order parameters merge correspondingly. There is region 
of temperatures in the graph, To < T < T* , where two different IRSB solutions 
coexist with the RS solution. The physical solution should have the largest free 
energy. Our calculations show that IRSB solution (with miRSB < 1) has larger 
free energy than the RS solution when Tq < T < Tirsb- So this IRSB solution 
is the physical one. The transition from the RS to the IRSB solution occurs 
jumpwise at Tirsb- At this point, _Frs = ^'irsb. So above Tirsb the RS-solution 
is the physical one- When T2 < T < Tirsb then the IRSB solution is also the 
physical one since it is stable with the respect to small perturbations having 2RSB 
symmetry- 



elements, we divide by mi/m2 groups of m2 elements each in order to find 2RSB 
solution. The parameter g"^ in 2RSB case we shall denote q™^^, if replicas a and 
P belong the same and at the same time the smallest subgroup (amount of gj"^^^ 
elements is equal to n(m2 — l)/2). The elements of q"^ will be labeled by q^'^^^ if 
replicase belong the same subgroup but this subgroup is not the smallest one (there 
are n{mi — m2)/2 replicas of that kind). Finally the notation ^q''^^ will be used if the 
replicas a and /? belong to different subgroups [there are n{n — mi)/2 replicas of this 
kind]- We set q'^^ = g"^ + Sq"^ and assume that the deviations of the IRSB 
solution from the 2RSB are small- We believe that one can neglect the changes of the 
order parameter WmsB ■ Then the free energy AF2 — i^2RSB — ^irsb can be conveniently 
represented as follows: 

AF2/NkT - ^ {{Sqrn'a + {i^^fh + [v^fc] , (33) 
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where the following substitution was used: 



Eq.(p3Dwas found in Ref. 27 for pairwise quadrupole interaction (J = 1, p = 
2)without any assumptions about the behavior of the order parameter Wmss- 

Parameters a, b, c, d can be expressed through averages found on IRSB-solution. 
For example, 

c = -2A(iKSB,,epi(m2 - mi)(l ~ mi)(l - maMp - l)(ri + vi)^P-^\ (36) 

Then we find that the parameter c < 0, if A(ip,sB,,„p, > 0, because from Eqs. (p^),(p^ 
follows that {ri + vi) > 0. Moreover the replica symmetry breaking scenario [0| 
assumes that in the limit n — > we obtain (7712 — mi) > 0, (1 — mi) > and 
(1 — 7712) > 0. Explicit form of the other parameters, namely, a,b,d, we give in the 
appendix (Appendix C). There it is shown that 5 < while A(iRs)rcpi > 0. 

It is shown in Appendix C that while U satisfies the reflection symmetry condition 
then the coefficient yi ^ and the stability is determined by the sign of A(ii,sB)rcpi 
(a < if A(iRSB)rcpi > 0). 

Direct numerical calculation in the subspace: U — quadrupole operator with 
J ~ 2, p ^ 3, shows that a < if A(ir,sB)rcpi > and at the same time (mi — 1) < 0, 
see Fig. I. Thus we can conclude that IRSB solution is stable indeed in this subgroup. 



3. Conclusions 



To conclude, we have demonstrated that IRSB solution in generalized p-spin models 
behaves differently depending on the symmetry of the operators U. If the operators U 
have zero trace for all odd powers, Tr t/'^^'^"'"^^ — for all integer k, then there trivial 
solution, = 0, for RS order parameter always exists. Therefore the bifurcation 
condition, A(RS)„p, — 0, can be satisfied for finite p > 3 only at Tq — 0. In this case 
IRSB solution smoothly branches only for p ~ 2. When p > 3 and Trt/(^'^+^) = 
IRSB solution can appear only discontinuously at the temperature defined by the 
condition, mi = 1 (if T^asB 7^ 0) as in usual p-spin model. IRSB solution is stable for 

-^(iRSB)ropl > 0- 

On the contrary when Tr {/(^fc+i) ^ Q solutions with broken replica symmetry 
may smoothly appear not only in the case of pair interactions, p = 2, but also at any 
finite p > 2. Point is that in this case the trivial RS-solution does not exist. The IRSB 
branching condition A(pis),„pi = is satisfied at finite temperature Tq. The solution 
with broken symmetry appears smoothly if mi(To) < 1 otherwise physically stable 
IRSB solution appears discontinuously at temperature where mi = 1. Algebraical 
expression for mi(To) is obtained. It depends only on RS-solution at Tq. 
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Appendix A. 

The only nonzero sums: 

lini (Sq"^f = -[r - (mi - l)vi]^ - mUl - miW; 



lim - y 5q°'Hq°'^ = [r - (rm - l)^!]^ 



(A.1) 
(A.2) 



lim i V ibq'^^f = -[r - (mi - VSv^f + 3mi(mi - l)[r - (mi - \)vyW + 



mi(mi - l)(2mi - l)?;i^; (A.3) 

lim - y 6q"^ 6q^^ Sq'^'^ = 2\r - (mi - l)vil'' - 3mi(mi - l)[r - (mi - l)viW - 



a,p,7 



mi^(mi - l)vi^; 



(A.4) 



lim - y (Sq^'^fSq"^ = [r - (mi - l)vi]^ - [r - (mi - l)t;i]t;i^ (A.5) 

a,/3,7 

hm - Sq^f^Sq^^Sq"^ = lim - 5q°'^h°'^h^^ = -[r - (mi - l)v{\^. (A.6) 



n->0 n 



n->0 n 



a,/3,7,(5 Q,/5,7,(5 

The prime on the sum means that only the superscripts belonging to the same 
5q are necessarily different in . 

Appendix B. 



w = 


'p{p- 
2 




L = 


p{p- 
2 


1)„(P-2)1^ 

HTiS 


C = 


p{p- 
2 


1)>-2)1' 


D = 


'p{p- 
2 


HRS 


and 






B2 = 


'p{p- 
2 


- l)„(p-2) 

'iRS 




'P{P- 
2 


" 1) (p-2) 

qns 



{{UfUi) - 2([>i'f/2f/3) + {U1U2U3U4)} : 
{{U^U2U3) - {U1U2U3U4)} ; 
{-{B2 + B'^) + 2B^ + B':,-Bi}- 
{-3B3 -B'^ + 3B4} ; 



(B.l) 
(B.2) 
(B.3) 
(B.4) 



6 ZD 
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B. = 



PjP - 1) (p-2) 



qks 



D 2 D 



-{(jlfj^UzUAUr,) 



t^B':^ = t^ 



PjP - 1) Jp-2) 



qks 



-{U1U2U3U4U5U6) + ^{U^U^UsUiU^) - \0iU2U3U4)- 

tVip-lf{p-2) (2p_5) 



{-{U!U2U3) + {U1U2USU4)} : 



p{p - 1) Jp-2) 



1 3 



^{UiU2U3UiU5U6) - {UlU2U3UiUr,) + l{ulu^u^Ui)+ 



1 ,,",3r>3\\ ^Mp-1)(P-2)jp-3) 



t^S4 = t^ 

+3{U^U^U3Ui) + 2{UlU2lJ3iJi) - 2{UIU2U3UaU5); 
K2 - {UfUi) - SiUfUiUs) + 21{UlUiU^U4) + 6{UfU2U3U4) - 

2{UlUlUl) + 10{UiU2U3U4U5U6} - 28{UlU2U3U4U5}; 

A2 = {U^U^)-{Uth 

A3 = {if!u2) - {U1U2U3); 

An = 0^) - mtill) + 2{UiUlUi): 

A12 = {utU2U3) + MuiulUs) - ifilui) - 'siuiulUsUi) - {u!u2); 
Ai3 = -{ufuiuD + {utuD + Hu!ulU3U4) - 



i~t^-w 



(B.5) 
(B.6) 
(B.7) 
(B.8) 
(B.9) 

2{U^UlU3) - {UtU2U3) - 20IU2U3U4U5) + 20IU2U3U4); (B.IO) 



Au = 3{UtUiU3U4) - 2{UtUiUz 



{UfU2U3) - 4([/2[/2C/3C/4C/5) + 4.{Ui U2U3U4) ■ 



(B.ll) 



iT.f \ 2*^P(P- 1) (p-2) 
= P ^ ^ <s 



+ 



P{P - 1) 



1 2 



p[r- (mi - l)vi] 



1 



,„(p-2)^(p-2)i - 
Wks Qrs + 

3P(P- l) >-2),2 . 
gas t A2 



p^[r - (mi - l)v{\ 



p{p - 1) 



p {[r - (mi - l)ui]^ + t;^mi(l - mi)} 



P{P- 1)]^ *''..,(p-2)>-3) 



Qrs A3 - 



48 



PjP - 1) (p-2)' 



p^[r- (mi - l)vi] — 



Ai 



p{p- 1) 



w'^r'\^-'^A,2 + 
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p {[r - (mi - l)v^f + v'rmil - rm)} ^ ~ 



p[r- (toi - l)vi]' 



r = 



It 



t^ p{p - 1) (p 



I pjp - 1) 
2 



where 



and 



Tr 



Tr[expS] 
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,„(p-2) 2(p-2) . 
Wrs fes ^13 + 



(B.12) 
(B.13) 

(B.14) 
(B.15) 



Appendix C. 

6 = -(1 - mi)2A(iRSB)repiP(p - l)(r-i + vi)^^-^^ - 



{l-m,r-p\p-lf{r,+v^f^^-'^ 



Ids' 



4 / dz^ 



-|mi Tr(;7^e''iKSB^ 

Tr e^i-RSB 



Tre''iBSi3 



Xr e^i^s-B 



-+ 



(mi - 4) y 



dz' 



Xr ffiiRSB 



Tr e*iRSB 



Xr e^iKs^ 



mi / d^; 





mi 




Tr e^^R^R 








Tre^iflS-B 



Xr e '^RSB 



2 \ 



> . 



(C.l) 



We shall further use the following notations: ^irsb = © and (Jq) = 
/dsG[Xree]'"\ 

It follows from the Cauchy-Schwarz inequality that the expression (Tri/^e®) (Tre®) > 

(iVf/e®)' follows from (E„ A.') (En^n") > (E„^ni3„)'. So 

^J^Jds^ (Tre®)"^"' (TrC/e«)' |(TrC/2e«) (Tre®) - (Tr;7e«)'| > 0. 
Then wc find similarly that / dxA{x) J dxB{x) > [/ dxA{x)B{x)Y, and 
m^j dz^'j^y ds"" {Tre^r^ I dsf (TVe®)"^"' (TVt^e®)'- 



y (is^ (Tre®)"^ (TVe®)"^"' (rrUe^y 



> 0. 



(C.2) 
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Then it follows that b < when A(iHSB)repi > 0. 
Here a = {ab — d?)/b, where 

a = -2mip{p - l)ri(f-2) } 1 - -p{p - l)n^P-^^ [{yi + - 2?/3) + 

4mit/4 - 2mi2/3 - imi^y^ - 4mi(l - TOi)y6 + mi(2 - mi)y2] } (C.3) 
d = -mi(l - mi)tVb - l)Vi(P-2)(ri + {2y7 - miyg - (2 - mOyg} , (C.4) 



and 
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